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Abstract. Quantum phase transitions occur at zero temperature when some non-thermal 
control-parameter like pressure or chemical composition is changed. They are driven by quan- 
tum rather than thermal fluctuations. In this review we first give a pedagogical introduction 
to quantum phase transitions and quantum critical behavior emphasizing similarities with 
and difi'erences to classical thermal phase transitions. We then illustrate the general concepts 
by discussing a few examples of quantum phase transitions occurring in electronic systems. 
The ferromagnetic transition of itinerant electrons shows a very rich behavior since the mag- 
netization couples to additional electronic soft modes which generates an effective long-range 
interaction between the spin fluctuations. We then consider the influence of rare regions on 
quantum phase transitions in systems with quenched disorder, taking the antiferromagnetic 
transitions of itinerant electrons as a primary example. Finally we discuss some aspects of 
the metal-insulator transition in the presence of quenched disorder and interactions. 
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1 Classical and quantum phase transitions 

1.1 Introduction 

Phase transitions have played, and continue to play, an essential role in shaping the 
world. The large scale structure of the universe is the result of a sequence of phase 
transitions during the very early stages of its development. Later, phase transitions 
accompanied the formation of galaxies, stars and planets. Even our everyday life is 
unimaginable without the never ending transformations of water between ice, liquid 
and vapor. Understanding phase transitions is a thus a prime endeavor of physics. 

Under normal conditions the phase transitions of water are so-called first-order 
transitions. They involve latent heat, i.e., a finite amount of heat is released while 
the material is cooled through an infinitesimally small temperature interval around the 
transition temperature. Phase transitions that do not involve latent heat, the so-called 
continuous transitions, are particularly interesting since the typical length and time 
scales of fluctuations of, e.g., the density, diverge when approaching the transition 
point. These divergences and the resulting singularities of physical observables are 
called the critical behavior. Understanding critical behavior has been a great challenge 
for theoretical physics. More than a century has gone by from the first discoveries until 
a consistent picture emerged. However, the theoretical concepts established during this 
development, viz., scaling and the renormalization group, now belong to the central 
paradigms of modern physics. 

The phase transitions we encounter in everyday life occur at finite temperature. 
These so-called thermal or classical^ phase transitions are driven by thermal fluctu- 
ations. In recent years a different class of phase transitions, the so-called quantum 
phase transitions, has started to attract a lot of attention. Quantum phase transitions 
occur at zero temperature when some non-thermal control parameter is changed. They 
are driven by quantum rather than thermal fluctuations. Quantum phase transitions 
in electronic systems have gained particular attention since some of the most exciting 



The justification for caUing all thermal phase transitions classical will become clear in Sec. |l.4| 
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discoveries in contemporary condensed matter physics, such as the locahzation prob- 
lem, various magnetic phenomena, integer and fractional quantum Hall effects, and 
high-temperature superconductivity are often attributed to quantum critical points. 

The purpose of this review is twofold. The first section gives a pedagogical in- 
troduction to the field of quantum phase transitions with a particular emphasis on 
the similarities with and the differences to classical thermal phase transitions. After 
brieffy sketching the historical development the basic concepts of continuous phase 
transitions and critical behavior are summarized. We then consider the question 'How 
important is quantum mechanics for the physics of phase transitions?' which leads 
directly to the distinction between classical thermal and quantum phase transitions. 
In the following sections these ideas are illustrated by discussing a number of exam- 
ples of quantum phase transitions occurring in electronic systems. Specifically, in Sec. 
a toy model for a quantum phase transition is considered, the so-called quantum 
spherical model. It can be solved exactly, providing an easily accessible example of 
a quantum phase transition. Sec. contains a discussion of the ferromagnetic quan- 
tum phase transition of itinerant electrons. It is demonstrated that the coupling of 
the magnetization to additional soft modes in the zero-temperature electron system 
changes the properties of the transition profoundly. The influence of disorder on quan- 
tum phase transitions is studied in Sec. ^ paying particular attention to rare disorder 
fluctuations. It is shown that they can change the universality class of the transition 
or even destroy the conventional critical behavior. In Sec. ^ we discuss some aspects of 
the metal-insulator transition of disordered interacting electrons. On the one hand we 
consider the influence of local moments on the transition by incorporating them into 
a transport theory. On the other hand we study the transition by means of large-scale 
numerical simulations. To do this, an efficient numerical method is developed, called 
the Hartree-Fock based diagonalization. It is shown that electron-electron interactions 
can lead to a considerable enhancement of transport in the strongly localized regime. 
Finally, Sec. g is devoted to a short summary and outlook. 

1.2 From critical opalescence to quantum criticality 

In 1869 Andrews 0] discovered a very special point in the phase diagram of carbon 
dioxide. At a temperature of about 31 °C and 73 atmospheres pressure the properties 
of the liquid and the vapor phases became indistinguishable. In the neighborhood 
of this point carbon dioxide strongly scattered light. Andrews called this point the 
critical point and the strong light scattering the critical opalescence. Four years later 
van der Waals ^J presented his doctoral thesis 'On the continuity of the liquid and 
gaseous states ' which contained one of the first theoretical explanations of critical 
phenomena based on the now famous van der Waals equation of state. It provides 
the prototype of a mean-field description of a phase transition by assuming that the 
individual interactions between the molecules are replaced by an interaction with a 
hypothetic global mean field. In the subsequent years similar behavior was found for 
many other materials. In particular, in 1895 Pierre Curie ||^ noticed that ferromagnetic 
iron also shows such a critical point which today is called the Curie point. It is located 
at zero magnetic field and a temperature of about 770 °C, the highest temperature for 
which a permanent magnetization can exist in zero field. At this temperature phases 
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differing by the direction of the magnetization become obviously indistinguishable. 
Again it was only a few years later when Weiss Q proposed the molecular-field theory 
of ferromagnetism which qualitatively explained the experiments. Like the van der 
Waals theory of the liquid-gas transition the molecular-field theory of ferromagnetism 
is based on the existence of a hypothetic molecular (mean) field. The so-called classic 
era of critical phenomena culminated in the Landau theory of phase transitions pj . 
Landau gave some very powerful and general arguments based on symmetry which 
suggested that mean-field theory is essentially exact. While we know today that this 
is not the case, Landau theory is still an invaluable starting point for the investigation 
of critical phenomena. 

The modern era of critical phenomena started when it was realized that there was 
a deep problem connected with the values of the critical exponents which describe how 
physical quantities vary close to the critical point. In 1945 Guggenheim g| realized that 
the coexistence curve of the gas-fiuid phase transition is not parabolic, as predicted by 
van der Waals' mean-field theory. At about the same time Onsager Q exactly solved 
the two-dimensional Ising model showing rigorously that in this system the critical 
behavior is different from the predictions of mean-field theory. After these observations 
it took about twenty years until a solution of the 'exponent puzzle' was approached. 
In 1965 Widom |8| formulated the scaling hypothesis according to which the singular 
part of the free energy is a generalized homogeneous function of the parameters. A 
year later, Kadanoff proposed a simple heuristic explanation of scaling based on the 
argument that at criticality the system essentially 'looks the same on all length scales'. 
The breakthrough came with a series of seminal papers by Wilson llfl] in 1971. He 
formalized Kadanoff's heuristic arguments and developed the renormalization group. 
For these discoveries, Wilson won the 1982 physics Nobel price. The development of 
the renormalization group initiated an avalanche of activity in the field which still 
continues. 

Today, thermal equilibrium phase transitions are well understood in principle, even 
if new interesting transitions, e.g., in soft condensed matter systems, continue to be 
found. In recent years the scientific interest has shifted towards new fields. One of these 
fields are phase transitions in non-equilibrium systems. They occur, e.g., in systems 
approaching equilibrium after a non-infinitesimal perturbation or in systems driven 
by external fields or non-thermal noise to a non-equilibrium (steady) state. Examples 
are provided by growing surfaces, chemical reaction-diffusion systems, or biological 
systems (see, e.g., Refs. [O-Oj). Non-equilibrium phase transitions are characterized 
by singularities in the stationary or dynamic properties of the non-equilibrium states 
rather than by thermodynamic singularities. 

Another very active avenue of research are quantum phase transitions which are the 
topic of this review. The investigation of quantum phase transitions was pioneered by 
Hertz [|lj] who built on earlier work by Suzuki |g5[ and Beal-Monod ijl^l . He developed 
a renormalization group method for magnetic transitions of itinerant electrons which 
was a direct generalization of Wilson's approach to classical transitions. He found that 
the ferromagnetic transition is mean-field like in all dimensions d > I. While Hertz' 
general scaling scenario at a quantum critical point is valid, his specific predictions for 
the ferromagnetic quantum phase transition are incorrect, as will be explained in Sec. 

I 
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In recent years quantum phase transitions in electronic systems have attracted 
considerable attention from theory as well from experiment. Among the transitions 
investigated in detail are the ferromagnetic transition of itinerant electrons, the anti- 
ferromagnetic transition associated with high-temperature superconductivity, various 
magnetic transitions in the heavy fermion compounds, metal-insulator transitions, 
superconductor-insulator transitions, and the plateau transition in quantum Hall sys- 
tems. This list is certainly incomplete and new transitions continue to be found. For 
reviews on some of these transitions see, e.g., Refs. [ pT-p2[ . There is also a very recent 
text book on quantum phase transitions by Sachdev [ 23[. 



1.3 Basic concepts of phase transitions and critical behavior 

Since the discoveries of scaling and the renormalization group a number of excellent 
text books on phase transitions and critical behavior have appeared (e.g., those by 
Ma pjl or Goldenfeld 123]). Therefore, in this section we only briefly collect the basic 
concepts which are necessary for the later discussion. 

A continuous phase transition can usually be characterized by an order parame- 
ter, a concept first introduced by Landau. An order parameter is a thermodynamic 
quantity that is zero in one phase (the disordered) and non-zero and non-unique in 
the other (the ordered) phase. Very often the choice of an order parameter for a 
particular transition is obvious as, e.g., for the ferromagnetic transition where the 
total magnetization is an order parameter. Sometimes, however, finding an appropri- 
ate order parameter is a complicated problem by itself, e.g., for the disorder-driven 
localization-delocalization transition of non-interacting electrons. 

While the thermodynamic average of the order parameter is zero in the disordered 
phase, its fluctuations are non-zero. If the phase transition point, i.e., the critical point, 
is approached the spatial correlations of the order parameter fluctuations become long- 
ranged. Close to the critical point their typical length scale, the correlation length ^, 
diverges as 

ecxi-'' (1) 

where v is the correlation length critical exponent and t is some dimensionless distance 
from the critical point. It can be defined as t = \T — Td/Tc if the transition occurs at 
a non-zero temperature Tc- In addition to the long-range correlations in space there 
are analogous long-range correlations of the order parameter fluctuations in time. The 
typical time scale for a decay of the fluctuations is the correlation (or equilibration) 
time Tc- As the critical point is approached the correlation time diverges as 

r, oc e" (X t-"' (2) 

where z is the dynamical critical exponent. Close to the critical point there is no 
characteristic length scale other than ^ and no characteristic time scale other than 
Tc!n As already noted by Kadanoff M, this is the physics behind Widom's scaling 
hypothesis, which we will now discuss. 



''Note that a microscopic xutofF scale must be present to explain non-trivial critical behavior, for 
details see, e.g., Goldenfeld [)25|. In a solid such a scale is, e.g., the lattice spacing. 
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Let us consider a classical system, characterized by its Hamiltonian 

H{pi,qi) = HkiniPi) + Hpot{qi) (3) 

where qi and pi are the generalized coordinates and momenta, and Hkin and Hpot 
are the kinetic and potential energies, respectively.^ In such a system 'statics and 
dynamics decouple', i.e., the momentum and position sums in the partition function 

^ = / n c^Ke-^-"/"-^ j W %e"^-'/'^'-^ = Zk^nZpot (4) 

are completely independent from each other. The kinetic contribution to the free 
energy density / = —{kBT/V) \ogZ will usually not display any singularities, since it 
derives from the product of simple Gaussian integrals. Therefore one can study the 
critical behavior using effective time-independent theories like the Landau-Ginzburg- 
Wilson theory. In this type of theories the free energy is expressed as a functional of 
the order parameter M(r) only. All other degrees of freedom have been integrated 
out in the derivation of the theory starting from a microscopic Hamiltonian. In its 
simplest form [^|l^, ^ valid, e.g., for an Ising ferromagnet, the Landau-Ginzburg- 
Wilson functional ^[M] reads 

$[Af] = / d'^r M{r) (-^r^+t\ M{r) + u j d'^r M^{r) - B j d'^r M{r), 
Z = /z?[Af]e-*[*^l , (5) 



where B is the field conjugate to the order parameter (the magnetic field in case of a 
ferromagnet). 

Since close to the critical point the correlation length is the only relevant length 
scale, the physical properties must be unchanged, if we rescale all lengths in the system 
by a common factor 6, and at the same time adjust the external parameters in such a 
way that the correlation length retains its old value. This gives rise to the homogeneity 
relation for the free energy density, 

f{t,B)^h-''f{th^'\Bby-). (6) 

Here ys is another critical exponent. The scale factor b is an arbitrary positive number. 
Analogous homogeneity relations for other thermodynamic quantities can be obtained 
by differentiating /. The homogeneity law (H) was first obtained phenomenologically 
by Widom |g]. Within the framework of the renormalization group theory it can be 
derived from first principles. 

In addition to the critical exponents i^, hb and z defined above, a number of other 
exponents is in common use. They describe the dependence of the order parameter 
and its correlations on the distance from the critical point and on the field conjugate 
to the order parameter. The definitions of the most commonly used critical exponents 
are summarized in Table ^ Note that not all the exponents defined in Table |l] are 
independent from each other. The four thermodynamic exponents a,/3, 7,(5 can all 

'^Velocity dependent potentials like in the case of charged particles in an electromagnetic field are 
excluded. 
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Table 1 Definitions of the commonly used critical exponents in the 'magnetic language', 
i.e., the order parameter is the magnetization m = (M), and the conjugate field is a magnetic 
field B. t denotes the distance from the critical point and d is the space dimensionality. (The 
exponent ys defined in (H) is related to S hy ys = d5/{l + 5).) 
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be obtained from the free energy (p) which contains only two independent exponents. 
They are therefore connected by the so-called scaling relations 



2- a = 2/? + 7 , 
2-a = I3{6 + 1) 



(7) 
(8) 



Analogously, the exponents of tlie correlation length and correlation function are con- 
nected by two so-called hyperscaling relations 



2 — a = dv ^ 
7 = (2-?7)i' 



(9) 
(10) 



Since statics and dynamics decouple in classical statistics the dynamical exponent z 
is completely independent from all the others. 

The critical behavior at a particular phase transition is completely characterized 
by the set of critical exponents. One of the most remarkable features of continuous 
phase transitions is universality, i.e., the fact that the critical exponents are the same 
for entire classes of phase transitions which may occur in very different physical sys- 
tems. These classes, the so-called universality classes, are determined only by the 
symmetries of the Hamiltonian and the spatial dimensionality of the system. This 
implies that the critical exponents of a phase transition occurring in nature can be 
determined exactly (at least in principle) by investigating any simplistic model system 
belonging to the same universality class, a fact that makes the field very attractive for 
theoretical physicists. The mechanism behind universality is again the divergence of 
the correlation length. Close to the critical point the system effectively averages over 
large volumes rendering the microscopic details of the Hamiltonian unimportant. 

The critical behavior at a particular transition is crucially determined by the rel- 
evance or irrelevance of order parameter fluctuations. It turns out that fluctuations 
become increasingly important if the spatial dimensionality of the system is reduced. 
Above a certain dimension, called the upper critical dimension d^, fluctuations are 
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irrelevant, and the critical behavior is identical to that predicted by mean- field the- 
ory (for systems with short-range interactions and a scalar or vector order parameter 
d^ = 4). Between dt^ and a second special dimension, called the lower critical di- 
mension d^, a phase transition still exists but the critical behavior is different from 
mean-field theory. Below d^ fluctuations become so strong that they completely sup- 
press the ordered phase. 

1.4 How important is quantum mechanics? 

The question of to what extent quantum mechanics is important for understanding a 
continuous phase transition is a multi-layered question. One may ask, e.g., whether 
quantum mechanics is necessary to explain the existence and the properties of the 
ordered phase. This question can only be decided on a case-by-case basis, and very 
often quantum mechanics is essential as, e.g., for superconductors. A different question 
to ask would be how important quantum mechanics is for the asymptotic behavior 
close to the critical point and thus for the determination of the universality class the 
transition belongs to. 

It turns out that the latter question has a remarkably clear and simple answer: 
Quantum mechanics does not play any role for the critical behavior if the transition 
occurs at a finite temperature. It does play a role, however, at zero temperature. In 
the following we will first give a simple argument explaining these facts. To do so it is 
useful to distinguish fiuctuations with predominantly thermal and quantum character 
depending on whether their thermal energy ksT is larger or smaller than the quantum 
energy scale HuJc- We have seen in the preceeding section that the typical time scale Tc 
of the fluctuations diverges as a continuous transition is approached. Correspondingly, 
the typical frequency scale ojc goes to zero and with it the typical energy scale 

hcUc OC \t\'''' . (11) 

Quantum fluctuations will be important as long as this typical energy scale is larger 
than the thermal energy ksT. If the transition occurs at some flnite temperature 
Tc quantum mechanics will thus become unimportant for |t| < t^ with the crossover 
distance tx given by 

tx OC T^/'^- . (12) 

We thus find that the critical behavior asymptotically close to the transition is en- 
tirely classical if the transition temperature Tc is nonzero. This justifies to call all 
finite-temperature phase transitions classical transitions, even if the properties of the 
ordered state are completely determined by quantum mechanics as is the case, e.g., 
for the superconducting phase transition of, say, mercury at Tc = 4.2 K. In these 
cases quantum fiuctuations are obviously important on microscopic scales, while clas- 
sical thermal fiuctuations dominate on the macroscopic scales that control the critical 
behavior. If, however, the transition occurs at zero temperature as a function of a 
non-thermal parameter like the pressure p, the crossover distance t^ — 0- (Note that 
at zero temperature the distance t from the critical point cannot be defined via the 
reduced temperature. Instead, one can define t = \p — Pc\/pc-) Thus, at zero tempera- 
ture the condition \t\ < t^ is never fulfilled, and quantum mechanics will be important 
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for the critical behavior. Consequently, transitions at zero temperature are called 
quantum phase transitions. 

Let us now generalize the homogeneity law (jm to the case of a quantum phase 
transition. We consider a system characterized by a Hamiltonian H. In a quantum 
problem kinetic and potential part of H in general do not commute. In contrast to 
the classical partition function (0) the quantum mechanical partition function does 
not factorize, i.e., 'statics and dynamics are always coupled'. The canonical density 
operator g-^/f'sT ^goks exactly like a time evolution operator in imaginary time r if 
one identifies 

l/keT = T = -iO/h (13) 

where Q denotes the real time. This naturally leads to the introduction of an imaginary 
time direction into the system. An order parameter field theory analogous to the 
classical Landau-Ginzburg- Wilson theory (H) therefore needs to be formulated in terms 
of space and time dependent fields. The simplest example of a quantum Landau- 
Ginzburg- Wilson functional, valid for, e.g., an Ising model in a transverse field, reads 



drjd'rM{v,r)[-—^-—^ 



$[Af] = ^ " dr j d'r M{v,T)\^-^-^+t]M{v,T) + 



+ u dT I d'^r M'^{y,t) - B I dr I d'^r M{t,t) . (14) 

Let us note that the coupling of statics and dynamics in quantum statistical dynamics 
also leads to the fact that the universality classes for quantum phase transitions are 
smaller than those for classical transitions. Systems which belong to the same classical 
universality class may display different quantum critical behavior, if their dynamics 
differ. 

The classical homogeneity law (p) for the free energy density can now easily be 
adopted to the case of a quantum phase transition. At zero temperature the imaginary 
time acts similarly to an additional spatial dimension since the extension of the system 
in this direction is infinite. According to (|2|), time scales like the zth power of a length. 
(In the simple example (nj) space and time enter the theory symmetrically leading to 
z = 1.) Therefore, the homogeneity law for the free energy density at zero temperature 
reads 

/(t, B) = h-^'^+'^f{t b^/'',B by'') . (15) 



Comparing (^_5|) and (|g) directly shows that a quantum phase transition in d dimen- 
sions is equivalent to a classical transition in d + z spatial dimensions. Thus, for a 
quantum phase transition the upper critical dimension, above which mean-field criti- 
cal behavior becomes exact, is reduced by z compared to the corresponding classical 
transition. 

Now the attentive reader may ask: Why are quantum phase transitions more than 
an academic problem? Any experiment is done at a non-zero temperature where, as 
we have explained above, the asymptotic critical behavior is classical. The answer is 
again provided by the crossover condition dl2): If the transition temperature Tc is 
very small quantum fluctuations will remain important down to very small i, i.e., very 



10 
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Fig. 1 Schematic phase diagram in the vicinity of a quantum critical point (QCP). The soUd 
Une marks the boundary between ordered and disordered phase. The dashed hues indicate 
the crossover between predominantly quantum or classical character of the fluctuations, and 
the shaded area denotes the critical region where the leading critical singularities can be 
observed. Paths (a) and (b) are discussed in the text. 



close to the transition. At a more technical level, the behavior at small but non-zero 
temperatures is determined by the crossover between two types of critical behavior, 
viz. quantum critical behavior at T = and classical critical behavior at non-zero 
temperatures. Since the 'extension of the system in imaginary time direction' is given 
by the inverse temperature l/ksT the corresponding crossover scaling is equivalent to 
finite size scaling in imaginary time direction. The crossover from quantum to classical 
behavior will occur when the correlation time Tc reaches l/fc^T which is equivalent to 
the condition (f2|). By adding the temperature as an explicit parameter and taking 
into account that in imaginary-time formalism it scales like an inverse time (Uq), we 
can generalize the quantum homogeneity law (HS) to finite temperatures. 



/(i, B, T) = b-^'^+^'^f{t b^/", B 6^« , T6" 



(16) 



The resulting phase diagram close to a quantum critical points will be of one of two 
qualitative different types. The first type describes situations where an ordered phase 
exists at finite temperature. These phase diagrams are illustrated in Fig. n. Here p 
stands for the (non-thermal) parameter which tunes the quantum phase transition. 
According to (12) the vicinity of the quantum critical point can be divided into re- 
gions with predominantly classical or quantum fluctuations. The boundary, marked 
by the dashed lines in Fig. |l|, is not sharp but rather a smooth crossover line. At 
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sufficiently low temperatures these crossover lines are inside the critical region (i.e., 
the region where the leading critical power laws can be observed). An experiment 
performed along path (a) will therefore observe a crossover from quantum critical be- 
havior away from the transition to classical critical behavior asymptotically close to 
it. At very low temperatures the classical region may become so narrow that it is 
actually unobservable in an experiment. 

In addition to the critical behavior at very low temperatures, the quantum critical 
point also controls the behavior in the so-called quantum critical region ||27| . This 
region is located at the critical p but, somewhat counter-intuitively, at comparatively 
high temperatures (where the character of the fluctuations is classical) . In this region 
the system 'looks critical' with respect to p but is driven away from criticality by 
the temperature (i.e., the critical singularities are exclusively protected by T). An 
experiment carried out along path (b) will therefore observe the temperature scaling 
at the quantum critical point. 

The second type of phase diagram occurs if an ordered phase exists at zero tem- 
perature only (as is the case for two-dimensional quantum antiferromagnets) . In this 
case there will be no true phase transition in any experiment. However, the system 
will display quantum critical behavior in the above-mentioned quantum critical region 
close to the critical p. 



2 Quantum spherical model 

2.1 Classical spherical model 

In the process of understanding a novel physical problem it is often very useful to 
consider a simple model which displays the phenomena in question in their most basic 
form. In the field of classical equilibrium critical phenomena such a model is the 
so-called classical spherical model which is one of the very few models in statistical 
physics that can be solved exactly but show non-trivial (i.e., non mean-field) critical 
behavior. The spherical model was conceived by Kac in 1947 in an attempt to simplify 
the Ising model. The basic idea was to replace the discrete Ising spins having only the 
two possible values Si = ±1 by continuous real variables between — oo and cx) so that 
the multiple sum in the partition function of the Ising model is replaced by a multiple 
integral which should be easier to perform. However, the multiple integral turned out 
to be not at all simple, and for a time it looked as if the spherical model was actually 
harder to solve than the corresponding Ising model. Eventually Berlin and Kac [g8| 
solved the spherical model by using the method of steepest descent to perform the 
integrals over the spin variables. Stanley [g9| showed that the spherical model, though 
created to be a simplification of the Ising model, is equivalent to the n — > oo limit of 
the classical n-vector modcl.^ Therefore, it can be used as the starting point for a 
1/n-expansion of the critical behavior. 

In the following years the classical spherical model was solved exactly not only for 



''in the classical n-vector model the dynamical variables are n-dimensional unit vectors. Thus, 
the Ising model is the 1-vector model, the classical XY-model is the 2-vector model and the classical 
Heisenberg model is the 3- vector model. 
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nearest neighbor ferromagnetic interactions but also for long-range power-law inter- 
actions [Q , random interactions |^ , 32 1 , systems in random magnetic fields ||3^, |4| , 
and disordered electronic systems with localized states [|5|. Moreover, the model has 
been used as a test case for the finite-size scaling hypothesis |^,|3^. Reviews on the 
classical spherical model were given by Joyce p8] and Khorunzhy et al. |39]. 

Because the classical spherical model possesses such a wide variety of applications 
in the field of classical critical phenomena, it seems natural to look for a quantum 
version of the model in order to obtain a toy model for quantum critical behavior. 
Actually, the history of quantum spherical models dates back at least as far as the 
history of quantum critical behavior. In 1972 Obermair Efl] suggested a canonical 
quantization scheme for a dynamical spherical model. However, this and later studies 
focused on the classical finite temperature critical behavior of the quantum model and 
did not consider the properties of the zero temperature quantum phase transition. 



2.2 Quantization of the spherical model 

The classical spherical model consists of N real variables Si G (— oo, oo) that interact 
with an external field hi and with each other via a pair potential Uij . The Hamiltonian 
is given by 

In order to make the model well-defined at low temperatures, i.e., in order to prevent 
a divergence of Si in the ordered phase, the values of Si are subject to an additional 
constraint, the so-called spherical constraint. Two versions of the constraint have been 
used in the literature, the strict and the mean constraints, defined by 

E^' = ^' (18) 

i 

J2iSf) - N, (19) 

i 

respectively. Here (. . .) is the thermodynamic average. Both constraints have been 
shown to give rise to the same thermodynamic behavior while other quantities like 
correlation functions differ. In the following we restrict ourselves to the mean spherical 
constraint which is easier to implement in the quantum case. The Hamiltonian ( |l7| ) 
has no internal dynamics. According to the factorization (^) it can be interpreted as 
being only the configurational part of a more complicated problem. Therefore, the 
construction of the quantum model consists of two steps: First we have to add an 
appropriate kinetic energy to the Hamiltonian which defines a dynamical spherical 
model which can be quantized in a second step. 

In order to construct the kinetic energy term we define canonically conjugate mo- 
mentum variables Pi which fulfill the Poisson bracket relations {Si, Pj} = Sij . The sim- 
plest choice of a kinetic energy term is then the one of Obermair Effl, Hkin — o^i-^i^ 
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where g can be interpreted as inverse mass. In this case, the complete Haniiltonian of 
the dynamical spherical model 

H = Hura + H,i = ^i{Ep^ + \T. ^'J^*^^ + E ^'^' + '^ ( E ^' - A (20) 

i i^j i \ i / 

is that of a system of coupled harmonic oscillators. Here we have also added a source 
term for the mean spherical constraint (E9|). (The value of fi has to be determined 
self-consistently so that (|ig| ) is fulfilled.) 

In order to quantize the dynamical spherical model (|20|) we use the usual canonical 
quantization scheme: The variables Si and Pt are reinterpreted as operators. The 
Poisson bracket relations are replaced by the corresponding canonical commutation 
relations 

[S„ Sj] = 0, [P,,Pj] = 0, and [S„ Pj] = ikS,, . (21) 

Equations ([l9|), (po|), and ( pi|) completely define the quantum spherical model. At 
large T or g the model is in its disordered phase (Si) = 0. The transition to an 
ordered state can be triggered by lowering g and/or T. 

It must be emphasized that this model does not mimic (or even describe) Heisen- 
berg-Dirac spins. Instead it is equivalent to the n -^ oo limit of a quantum rotor 
model which can be seen as a generalization of an Ising model in a transverse field. Of 
course, the choices of the kinetic energy and quantization scheme are not unique. In 



agreement with the general discussion in Sec. 1.4 different choices will lead to different 
critical behavior at the quantum phase transition, while the classical critical behavior 
is the same for all these models. An example of a different quantization of the spherical 
model was given by Nieuwenhuizcn [Q. It leads to a dynamical behavior that more 
closely resembles that of Heisenberg-Dirac spins than our choice. For a more detailed 
discussion of these questions see also Ref . Q . 

2.3 Quantum phase transitions 

The quantum spherical model defined in eqs. (|l^), (|o|), and (|2|) can be solved exactly 
since it is equivalent to a system of coupled harmonic oscillators. This was done in 
Ref. [Q for a model with arbitrary translationally invariant interactions (long-range 
as well as short-range) in a spatially homogeneous external field. The resulting free 
energy reads 

ksT h^ kgT v^ / u;(k) \ , , 

/ = ^—\nZ = -n + -2— V n 2sinh-r^-- , (22) 

with w(fc) given by cu^ = 2g[^ + U(k)/2], where U{k) is the Fourier transform of the 
interaction Uij . The spherical constraint which determines fj, is given by 

„ 9f /i^ 1 v-^ q , uj(k) , , 

«=ai^-i+V + ivE4:)-th^. (23) 

As usual in spherical models the critical behavior is determined by the properties of 
the solutions of (E3) for small /i. At any finite temperature the coth-term can be 
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Table 2 Critical exponents at the quantum and classical phase transitions of the quantum 
spherical model as functions of the dimensionality d and the exponent x which character- 
izes the long-wavelength behavior of the interaction U{k) ~ \k\^ (short-range interactions 
correspond to a:: = 2). 





Quantum transition 
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expanded giving the same leading long-wavelength and low frequency terms as in the 
classical spherical model (p7|). As expected, the resulting critical behavior at finite 
temperatures is therefore that of the classical spherical model. 

At zero temperature, the coth-tcrm in (|2^) is identical to one. Thus, the leading 
long- wavelength and low frequency terms are different from the classical case. This 
gives rise to the quantum critical behavior being different from the classical one. If 
the interaction Uij in the Hamiltonian is short ranged, the dynamical exponent turns 
out to be z = 1. For a power-law interaction, parameterized by the singularity of the 
Fourier transform of the interaction, Uk c>c |fc|^ for /c — > 0, we obtain z = x/2. In both 
cases the quantum critical behavior of the d-dimensional quantum spherical model is 
the same as the classical critical behavior of a corresponding d + z-dimensional model. 
The critical exponents for the quantum and classical phase transitions are summarized 
in Table |. 

In order to describe the crossover between the quantum and classical critical be- 
haviors the crossover scaling form of the equation of state was derived. This is only 
possible below the upper critical dimension. Above, crossover scaling breaks down. 
This is analogous to the breakdown of finite-size scaling in the spherical model above 
the upper critical dimension. It can be attributed to a dangerous irrelevant variable. 

In Ref. p3| the influence of a quenched random field on the quantum phase transi- 
tion was considered. The quantum spherical model can be solved exactly even in the 
presence of a random field without the necessity to use the replica trick. It was found 
that the quantum critical behavior is dominated by the static random field fluctua- 
tions rather than by the quantum fluctuations. Since the random field fluctuations are 
identical at zero and finite temperatures it follows that in the presence of a random 
field quantum and classical critical behavior are identical. 
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3 Ferromagnetic quantum phase transition of itinerant electrons 

3.1 Itinerant ferroniagnets 

In the normal metallic state the electrons form a Fermi liquid, a concept introduced 



by Landau 1 44 , 45| . In this state the excitation spectrum is very similar to that of 
a non-interacting Fermi gas. The basic excitations are weakly interacting fermionic 
quasiparticles which behave like normal electrons but have renormalized parameters 
like an effective mass. However, at low temperatures the Fermi liquid is potentially 
unstable against sufficiently strong interactions, and some type of a symmetry-broken 
state may form. This low-temperature phase may be a superconductor, a charge 
density wave, or a magnetic phase, e.g., a ferromagnet, an anti-ferromagnet, or a 
spin glass, to name a few possibilities. In general, it will depend on the microscopic 
parameters of the material under consideration what the nature of the low-temperature 
phase and, specifically, of the ground state will be. Upon changing these microscopic 
parameters at zero temperature, e.g., by applying pressure or an external field or by 
changing chemical composition, the nature of the ground state may change, i.e., the 
system may undergo a quantum phase transition. 

In this Section we will discuss a particular example of such a quantum phase 
transition, viz. the ferromagnetic quantum phase transition of itinerant electrons. 
Most of the Section will be devoted to clean itinerant electrons but we will also briefly 
consider the influence of disorder on the ferromagnetic transition. 

The experimentally best studied example of a ferromagnetic quantum phase tran- 
sition of itinerant electrons is probably provided by the pressure-tuned transition in 
MnSi M,M . MnSi belongs to the class of so-called nearly or weakly ferromagnetic ma- 
terials. This group of metals, consisting of transition metals and their compounds such 
as ZrZn2, TiBe2, NisAl, and YC02 m addition to MnSi are characterized by strongly 
enhanced spin fluctuations. Thus, their ground state is close to a ferromagnetic in- 
stability which makes them good candidates for actually reaching the ferromagnetic 
quantum phase transition in experiment by changing the chemical composition or 
applying pressure. 

At ambient pressure MnSi is paramagnetic for temperatures larger than Tc — 30 K. 
Below Tc it orders magnetically. The order is, however, not exactly ferromagnetic but 
a long- wavelength (190 A) helical spin spiral along the (111) direction of the crystal. 
The ordering wavelength depends only weakly on the temperature, but a homogeneous 
magnetic field of about 0.6 T suppresses the spiral and leads to ferromagnetic order. 
One of the most remarkable findings about the magnetic phase transition in MnSi is 
that it changes from continuous to first order with decreasing temperature as is shown 
in Fig. ^. Specifically, in an experiment carried out at low pressure (corresponding to 
a comparatively high transition temperature) the susceptibility shows a pronounced 
maximum at the transition, reminiscent of the singularity expected from a continu- 
ous phase transition. In contrast, in an experiment at a pressure very close to (but 
still smaller than) the critical pressure the susceptibility does not show any sign of a 
divergence at the phase transition. Instead, it displays a finite step suggestive of a 
first-order phase transition. 

A related set of experiments is devoted to a phenomenon called the itinerant elec- 
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iron metamagnetism. Here a high magnetic field is applied to a nearly ferromagnetic 
material such as Co(Sei_a;Sa;)2 iQ or Y(Coi_a;Al2:)2 |Q. At a certain field strength 
the magnetization of the sample shows a pronounced jump. This can easily be ex- 
plained if we assume that the free energy as a function of the magnetization has the 
triple-well structure characteristic of the vicinity of a first-order phase transition. In 
zero field the side minima must have a larger free energy than the center minimum 
(since the material is paramagnetic in zero field). The magnetic field essentially just 
"tilts" the free energy function. If one of the side minima becomes lower than the 
center (paramagnetic) one, the magnetization jumps. 

In the literature the first-order transition in MnSi at low temperatures as well 
as the itinerant electron metamagnetism have been attributed to sharp structures in 
the electronic density of states close to the Fermi energy which stem from the band 
structure of the particular material. These structures in the density of states can 
lead to a negative quartic coefficient in a magnetic Landau theory and thus to the 
above mentioned triple-well structure. In the next section it will be shown, however, 
that the two phenomena are generic since they are rooted in the universal many-body 
physics underlying the transition. Therefore, they will occur for all nearly or weakly 
ferromagnetic materials irrespective of special structures in the density of states. 
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3.2 Landau-Ginzburg-Wilson theory of the ferromagnetic quantum phase transition 

From a theoretical point of view, the ferromagnetic transition of itinerant electrons 
is one of the most obvious quantum phase transitions. It was also one of the first 



quantum phase transitions investigated in some detail. Hertz 14 1 studied a simple 
microscopic model of interacting electrons and derived a Landau-Ginzburg-Wilson 
theory for the ferromagnetic quantum phase transition. Hertz then analyzed this the- 
ory by means of renormalization group methods which were a direct generalization of 
Wilson's treatment of classical transitions. He found a dynamical exponent of 2; = 3. 



According to the discussion in Sec. 1.4 this effectively increases the dimensionality of 
the system from d to d + 3. Therefore, the upper critical dimension of the quantum 
phase transition would be (i+ — 1, and Hertz concluded that the critical behavior 
of the ferromagnetic quantum phase transition is mean-field like in all physical di- 
mensions d > 1. While it was later found [Sy] that Hertz' description of the finite 
temperature phenomena close to the quantum critical point was incomplete, it was 
generally believed that the main qualitative results of his model at zero temperatures 
apply to real itinerant ferromagnets as well.^ However, in 1994 Sachdev [£2[ showed 
that Hertz' results in dimensions below one (an academic but still interesting case) 
violate an exact exponent equality. 

Vojta, Belitz, Narayanan, and Kirkpatrick |53] have revisited the ferromagnetic 
transition of itinerant electrons. They have shown that the properties of the transition 
are much more complicated since the magnetization couples to additional, non-critical 
soft modes in the electronic system. Mathematically, this renders the conventional 
Landau-Ginzburg-Wilson approach invalid since an expansion of the free energy in 
powers of the order parameter does not exist. Physically, the additional soft modes 
lead to an effective long-range interaction between the order parameter fluctuations. 
This long-range interaction, in turn, can change the character of the transition from a 
continuous transition with mean-field exponents to either a continuous transition with 
non-trivial (non mean-field) critical behavior or even to a first order transition. 

The derivation of the order parameter field theory [5^,0 follows Hertz ]1J] in spirit, 
but the technical details are considerably different. Let us consider a microscopic 
model Hamiltonian H = Hq + iJox of interacting fermions. iJox is the exchange 
interaction which is responsible for the ferromagnetism, iJg does not only contain the 
free electron part but also all interactions except for the exchange interaction. Using 
standard manipulations (see, e.g., Ref. |5^) the partition function is written in terms 
of a functional integral over fermionic (Grassmann) variables. After introducing the 
magnetization field M(r,r) via a Hubbard-Stratonovich transformation [p6|,p7| of the 
exchange interaction, a cumulant expansion is used to integrate out the fermionic 
degrees of freedom. The partition function Z takes the form 



_ ^-Fo/T 



D[M\ exp[-$[M]] , (24) 



''In order to obtain a quantitative description Moriya and Kawabata developed, a more sophisticated 
theory, the so-called self-consistent renormalization theory of spin fluctuations felj] . This theory is very 
successful in describing magnetic materials with strong spin fluctuations outside the critical region. Its 
results for the critical behavior at the ferromagnetic quantum phase transition are, however, identical 
to those of Hertz. 
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where Fq is the non-critical part of the free energy. With the four-vector notation 
with X — (x, r) and J dx = J dx L dr the resulting Landau-Ginzburg- Wilson 

free energy functional reads 



$[M] = - dxdy 



—5{x-y)-x^'^\x~y) 
1 1 



M(x) • M(y) -f (25) 



n=3 "^' •' 

where Ft is the spin-triplet (exchange) interaction strength. The coefRcients in the 
Landau-Ginzburg- Wilson functional are the connected n-point spin density correlation 
functions x'"' of the reference system Hq which is a conventional Fermi liquid. The 
famous Stoner criterion [ p8[ of ferroniagnetisni, Ft g(eF) > 1 (here giep) is the density 
of states at the Fermi energy) can be rediscovered from the stability condition of the 
Gaussian term of <&[M], if one takes the spin susceptibility x^^'' to be that of non- 
interacting electrons (in which case x'^^H'i ^ 0, il = 0) = g{eF))- 

The long- wavelength and long-time properties of the spin-density correlation func- 
tions of a Fermi liquid were studied [p9| using diagrammatical perturbation theory 
in the interaction. Somewhat surprisingly, all these correlation functions generically 
(i.e., away from any critical point) show long-range correlations in real space which 
correspond to singularities in momentum space in the long- wavelength limit q ^ 0. 
While analogous generic long-range correlations in time (the so-called long-time tails) 
are well known from several interacting systems, long-range spatial correlations in 
classical systems are impossible due to the fluctuation-dissipation theorem. They are 
known, however, in non-equilibrium steady states (see, e.g., Ref. [|60[ ). The physi- 
cal reason for the singularities in the coefRcients x*^"-* of the Landau-Ginzburg- Wilson 
functional is that in the process of integrating out the fermionic degrees of freedom 
the soft particle-hole excitations have been integrated out, too. It is well known from 
classical dynamical critical phenomena ||61| that integrating out soft modes leads to 
singularities in the resulting effective theory. 

Specifically, it was found [B9| that the static spin susceptibility x^^H^") behaves like 
j,-(2d-i) jTqj. igji-gQ distances r. The leading long- wavelength dependence therefore has 
the form 

X(2)(q)/X(^)(0) = l + Crf(|q|/2fcf)'^-i + 0(|q|2) (d < 3) (26) 

while in d = 3 the non-analyticity takes the form C3(|q|/2fci?)^ ln(2fci?/|q|). Here kp 
is the Fermi momentum and Cd and C3 are dimensionless constants. Note that these 
singularities only exist at zero temperature and in zero magnetic field since both a 
finite temperature and a magnetic field give the particle-hole excitations a mass. 
Using (^6|), and with / = ^ ^^X^iiji the Gaussian part of $ can be written, 

^^^^[M]= fM{q)[to + Cd\cif+C2q' + cn\n\/\q\]Mi-q) . (27) 

Jq 

Here to = 1 — Ftx'-^Hl ~^ 0,w„ = 0) is the bare distance from the critical point, and 
Co is another constant. Physically, the non- analytic term in the Gaussian part of <f> 
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represents a long-range interaction of the spin fluctuations which is self-generated by 
the electronic system. For the same physical reasons for which the non-analyticity 
occurs in x^^^ the higher coefficients x^""* (" > 2) in (25) in general diverge for 



zero frequencies and wave numbers. Consequently, the free energy functional (Gl 
is mathematically ill-defined. However, it will nonetheless be possible to extract a 
considerable amount of information. 

The sign of the non-analyticity in the Gaussian term merits some attention since 
it will be responsible for the qualitative features of the ferromagnetic quantum phase 
transition. Perturbation theory to second order in Tt yields Cd < [p9|. This is the 
generic case, and it is consistent with the well-known notion that correlation effects in 
general decrease the effective Stoner coupling |6^]. However, Ref. [^ has given some 
possible mechanisms for Cd to be positive at least in some materials. 

3.3 Phase transition scenarios 

Depending on the sign of the non-analyticity in the Gaussian term ( |27| ) of the free 
energy functional there will be different scenarios for the ferromagnetic quantum phase 
transition |p3]. 

We first discuss the generic case of c^ < 0. Here the free energy reduces with 
increasing q from zero which implies that a continuous transition to a ferromagnetic 
state is impossible at zero temperature. Two possible scenarios for the phase transition 
arise for Cd < 0- The first scenario is based on the observation that a finite thermody- 
namic magnetization m = (|M(a;)|), which acts similarly to a magnetic field, cuts off 
the singularities in the coefficients of the order parameter field theory. Therefore, the 
non-analyticity in x leads to an analogous non-analyticity in the magnetic equation 
of state, which takes the form 

tm ~ Vdm'^ + um^ = H (rf < 3) , (28) 

tm-V3m^ln{l/m) + um^ = H (d = 3) , (29) 

where t tunes the transition and u, Vd and V3 are positive constants. H denotes the 
external magnetic field. This equation of state describes a first-order phase transition 
since the next-to-leading term for small m has a negative sign. This scenario was 
investigated in some detail in Ref. |Q]. Since the non-analyticities in x'^' and the 
equation of state are cut off by a finite temperature, the transition will be of first 
order at very low T but turn second order at higher temperatures. The two regimes 
are separated by a tricritical point. This is exactly the behavior found experimentally 
inMnSi ^^. 

The second possible scenario for the quantum phase transition arising if c^; < 
is that the ground state of the system will not be ferromagnetic but instead a spin- 
density wave at finite q. This scenario has not been studied in much detail so far, 
but work is in progress. It is tempting to interpret the spiral ordering in MnSi as a 
signature of this finite-g instability. This is, however, not very likely since a finite-g 
instability caused by the long-range interaction will be strongly temperature dependent 
due to the temperature cutoff of the singularities. As mentioned above, experimentally 
the ordering wave vector is essentially temperature independent. Further work will be 
necessary to decide which of the two possible scenarios, viz. a first-order ferromagnetic 
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transition or a continuous transition to modulated magnetic order, is realized under 
what conditions. Moreover, let us point out, that in d = 3 the non-analyticity is only 
a logarithmic correction and would hence manifest itself only as a phase transition at 
exponentially small temperatures, and exponentially large length scales. Thus, it may 
well be unobservable experimentally for some materials. 

We now turn to the second case, Cd > which can happen, if one of the conditions 
discussed in Rcf. M] is fulfilled. In this case the self-generated long-range interaction 
is ferromagnetic. Consequently, the ferromagnetic quantum phase transition will be 
a conventional second order phase transition, which can be analyzed by standard 
renormalization group methods. A tree level analysis shows that the Gaussian theory 
is sufficient for dimensions d > d^ = 1 since all higher order terms are irrelevant. 
We are therefore able to obtain the critical behavior exactly, yet due to the long-range 
interaction it is not mean field-like. The results of this analysis |54[| can be summarized 
as follows. At zero temperature the equation of state close to the quantum critical 
point reads 

tm + Vdm!^ -f um^ ^ H (rf < 3) , (30) 

tm + vj,m^ \n{\ / m) + urrC' = H (rf = 3) , (31) 

Again, u and v are positive constants. Note the different sign of the non-analytic 
terms compared to eqs. (Eq, ES). From the equation of state one obtains the criti- 
cal exponents (3 and 5 while the correlation length exponent v, the order parameter 
susceptibility exponent 77, and the dynamical exponent z can be directly read of the 
Gaussian part of $, eq. (pTj). We find j3 = v = l/{d — 1), 'r] = 3 — d, S~z — d 
for 1 < d < 3. These exponents 'lock into' mean-field values fi ^ v = 1/2, r/ = 0, 
5 — z = i for d > 3. In d = 3, there are logarithmic corrections to power-law scaling. 
At finite temperature, we find homogeneity laws for ni, and for the magnetic sus- 
ceptibility, Xm, 

m(t, T, H) = b-f^/^mith^/", Tb'^^", Hb^^^'") , (32) 

X™(t,T,i7) = 6^/''x™(i6^/^T6^/^mW-) , (33) 

where b is an arbitrary scale factor. The susceptibility exponent 7 and the crossover 
exponent (j) that describes the crossover from the quantum to the classical Heisenberg 
fixed point (FP) are given by 7 = /3((5 — 1) = 1, </) = i^ for all d > 1. Notice that the 
temperature dependence of the magnetization is not given by the dynamical exponent. 
However, z controls the temperature dependence of the specific heat coefficient, 71/ = 
cy jT , which has a scale dimension of zero for all d, and logarithmic corrections to 
scaling for all d < 3 

7v(t,r,i/)==e(3-d) ln6-K7v(t6l/^T6^mW'^) . (34) 

The singularities in the spin density correlation functions do not only influence the 
properties of the quantum phase transition but also those of the ferromagnetic phase. 
An example is the dispersion relation of the ferromagnetic spin waves [p5[ . Since the 
non-analyticities are cut off by a finite magnetization it turns out that the disper- 
sion relation remains w oc (7^ but the prefactor picks up a non-trivial magnetization 
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Fig. 3 Phase diagrams of disordered itinerant ferromagnets in the T-t-plane showing first 
order (solid lines) and second order (dashed lines) transitions. G is a dimensionless measure 
of the disorder, and CP, CEP, TCP, and TCEP denote critical points, critical endpoints, 
tricritical points and tricritical endpoints, respectively. The first panel corresponds to the 
experimental results on MnSi M,[47|. See Ref. [b4| for more details. 



dependence (different from being proportional to m as in Stoner theory). For small 
magnetization m we find 
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(35) 
(36) 



Until now, the corrections to mean-field theory predicted by (p3) and (pm have not 
been observed experimentally. 



3.4 Influence of disorder 

In this subsection we briefly discuss the influence of quenched non-magnetic disorder 
on the ferromagnetic quantum phase transition. An approach similar to that of the 
clean case has also been developed for the dirty case BO], and the resulting effective 
theory is very similar. Again, the magnetization couples to additional soft modes 
(here with diffusive dynamics) which leads to an effective long-range interaction. The 
singularities are even stronger than in the clean case, but they have the opposite sign 
so that the self-generated long-range interaction is generically ferromagnetic. Thus in 
the presence of disorder there will be a competition between the ballistic and diffusive 
singularities, and the temperature which cuts off both. For weak disorder the first- 
order transition will survive, while larger disorder leads to a continuous transition. As 
shown in Ref. |64[| , the phase diagram becomes very rich, showing several multicritical 
points and even regions with metamagnetic behavior (see Fig. |^). The properties of 
the continuous quantum phase transition occurring for stronger disorder can again be 
analyzed by standard renormalization group methods. It turns out that as in the clean 
case the Gaussian theory is sufficient since all higher order terms are irrelevant. The 
resulting critical exponents are 7 = 1 for all d > 2 and ly = l/{d — 2), ry = 4 — d, z ~ d 
for 2 < d < 4. These exponents lock into their mean field values i^ = 1/2, 77 = 0, and 
z = 4 for d > 4. In addition to d = 4, d = 6 also plays the role of an upper critical 
dimension, and one has (3 = 2/((i — 2), (5 = d/2 for 2 < d < 6, while (3 = 1/2, (5 = 3 
for d > 6. 
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4 Influence of rare regions on magnetic quantum phase transitions 

4.1 Disorder, rare regions, and the GrifEths region 

The influence of static or quenched disorder on the critical properties of a system near 
a continuous phase transition is a very interesting problem in statistical mechanics. 
While it was initially suspected that quenched disorder always destroys any critical 
point pj, this was soon found to not necessarily be the case. Harris |6^, ^ found 
a convenient criterion for the stability of a given critical behavior with respect to 
quenched disorder: If the correlation length exponent i' obeys the inequality i' > 2/d, 
with d the spatial dimensionality of the system, then the critical behavior is unaffected 
by the disorder. In the opposite case, i' < 2/d, the disorder modifies the critical 
behavior | |70| . This modification may either (i) lead to a new critical point that has a 
correlation length exponent v > 2/d and is thus stable, or (ii) lead to an unconventional 
critical point where the usual classification in terms of power-law critical exponents 
looses its meaning, or (iii) lead to the destruction of a sharp phase transition. The first 
possibility is realized in the conventional theory of random- Tc classical ferromagnets 
[ 37 1 , and the second one is probably realized in classical ferromagnets in a random field 
[ 71 - |73| . The third one has occasionally been attributed to the exactly solved McCoy- 
Wu model |74|-|7q| . This is misleading, however, as has recently been emphasized in 
Ref. |77|| ; there is a sharp, albeit unorthodox, transition in that model, and it thus 
belongs to category (ii). 

Independent of the question of if and how the critical behavior is affected, disorder 
leads to very interesting phenomena as a phase transition is approached. Disorder 
in general decreases the critical temperature Tc from its clean value T^. In the tem- 
perature region Tc < T < T^ the system does not display global order, but in an 
infinite system one will find arbitrarily large regions that are devoid of impurities, 
and hence show local order, with a small but non-zero probability that usually de- 
creases exponentially with the size of the region. These static disorder fluctuations 
are known as 'rare regions', and the order parameter fluctuations induced by them 
as 'local moments' or 'instantons'. Since they are weakly coupled, and flipping them 
requires to change the order parameter in a whole region, the local moments have very 
slow dynamics. Griffiths JTg] was the first to show that they lead to a non-analytic free 
energy everywhere in the region Tc < T < T^ , which is known as the Griffiths phase, 
or, more appropriately, the Griffiths region. In generic classical systems this is a weak 
effect, since the singularity in the free energy is only an essential one. An important 
exception is the McCoy- Wu model [£4| , which is a 2D Ising model with bonds that are 
random along one direction, but identical along the second direction. The resulting 
infinite-range correlation of the disorder in one direction leads to very strong effects. 
As the temperature is lowered through the Griffiths region, the local moments cause 
the divergence of an increasing number of higher order susceptibilities, d^m/dB"' 
(n > 2), starting with large n. Even the average susceptibility proper, x^^^ ^ dm/dB, 
diverges at a temperature T,^ > Tc, although the average order parameter does not 
become non-zero until the temperature reaches Tc. This is caused by rare fiuctuations 
in the susceptibility distribution, which dominate the average susceptibility and make 
it very different from the typical or most probable one. 
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Surprisingly little is known about the influence of the Grifliths region and related 
phenomena on the critical behavior. Recent work 179] on a random- Tc classical Ising 
model has suggested that it can be profound, even in this simple model where the 
conventional theory predicts standard power-law critical behavior, albeit with critical 
exponents that are different from the clean case. The authors of Ref. [f9[ have shown 
that the conventional theory is unstable with respect to perturbations that break the 
replica symmetry. By approximately taking into account the rare regions, which are 
neglected in the conventional theory, they found a new term in the action that actually 
induces such perturbations. In some systems replica symmetry breaking is believed to 
be associated with activated, i.e. non-power law, critical behavior. Reference 179] thus 
raised the interesting possibility that, as a result of rare-region effects, the random- Tc 
classical Ising model shows activated critical behavior, as is believed to be the case for 
the random-field classical Ising model 1 71 - 7^] , although in the case of the random- Tc 



model no final conclusion about the fate of the transition could be reached. 

Griffiths regions also occur in the case of quantum phase transitions (for an ex- 
perimental example see Ref. [pQ]). Their consequences for the critical behavior are 
even less well investigated than in the classical case, with the remarkable exception 
of certain ID systems. Fisher |77[ ] has investigated quantum Ising spin chains in a 
transverse random field. These systems are closely related to the classical McCoy- Wu 
model, with time in the quantum case playing the role of the 'ordered direction' in 
the latter. He has found activated critical behavior due to rare regions. This has 
been confirmed by numerical simulations ] pl[ . Other recent simulations ]p2| suggest 
that this type of behavior may not be restricted to ID systems, raising the possibility 
that exotic critical behavior dominated by rare regions may be generic in quenched 
disordered quantum systems, independent of the dimensionality and possibly also of 
the type of disorder. 



4.2 Itinerant quantum antiferromagnets 

Within the conventional theory m% of critical behavior in systems with quenched 
disorder the first step consists of averaging over the disorder, usually via the replica 
trick ]Q. The resulting effective theory is then analyzed perturbatively. However, 
the rare regions are a non-perturbative effect since the probability for their occurrence 
is exponentially small in the disorder strength. Therefore, rare regions are neglected 
within the conventional theory. 

Narayanan, Vojta, Belitz, and Kirkpatrick |8J] have developed a generalization 
of the conventional theory of quantum phase transitions in the presence of quenched 
disorder. This theory, which is similar to that of Ref. m^ for classical transitions, 
includes the effects of the rare regions. The basic idea is not to average over the 
disorder at the beginning but to work with a particular disorder configuration until 
the rare regions are identified. Only after their effects have been incorporated into the 
theory, the disorder average is carried out. 

In the following we illustrate this theory taking the itinerant quantum antifer- 
romagnet as the primary example. The starting point is the order parameter field 
theory for the itinerant quantum antiferromagnet derived by Hertz |L4| . The Landau- 
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Ginzburg- Wilson free energy functional reads 

S* = dx dy (/)(x) ro(x, y) 4>{y) + u dx ((/)(a;) • (j>ix)f , (37) 



where (p is the staggered magnetization. T(x,y) is the bare two-point vertex function, 
whose Fourier transform is 

ro(q,l^„) = (to + q2 + |w„|) . (38) 

Disorder is introduced by making the distance t from the critical point a random 
function of position, i(x) = io + ^t{'x.), where dt{x.) obeys a Gaussian distribution with 
zero mean and variance A. 

Instead of averaging over the disorder we now determine saddle point solutions of 



the unaveraged Landau-Ginzburg- Wilson functional (37). Due to the disorder there 
will be spatial regions in which the system wants to order (t(x) < 0) even if it is globally 
in its disordered phase (io > 0). These rare regions or islands will support locally 
nonzero saddle-point solutions. Outside of the islands, the solution is exponentially 
small. Thus, the islands are effectively decoupled. For a system with N islands, and in 
the case of Ising symmetry, there will be 2^ almost degenerate saddle-point solutions 
that can be constructed by considering all possible distributions of the sign of the 
order parameter on the islands. For a continuous order parameter symmetry there is 
a whole manifold of almost degenerate saddle points. This complicated structure of 
the free energy landscape is responsible for the failure of the conventional theory as is 
known from the random field Ising model [|7l|-[73|. 

Now, the crucial point is that for a complete theory one has to take into account 
fluctuations around all of these saddle points. As was shown in Ref. [Q the saddle 
point configurations act as an additional source of disorder in the system. Since 
the saddle points are time-independent this disorder is static, but it is self-generated 
and thus in equilibrium with the rest of the system. Therefore, taking into account 
fluctuations around all saddle points leads to the appearance of static annealed disorder 
in addition to the underlying quenched disorder. (Some general aspects concerning 
annealed disorder and quantum phase transitions are discussed in Ref. ^^■) At this 
point in the calculation the average of the quenched disorder is carried out by means 
of the replica trick. The resulting effective theory for the fluctuations ip takes the form 

5eff = ^ /'dxdyv'"(a^)ro(x,y)¥P"(y)+7i^ /da; (¥P"(a;).(p"(a;))2 

- Aj2fdxdyS{^-y){^^{x)f{ip^\y))' 

- TwY,fdxdyS{K~y){^^{x))'{^"{y)f. (39) 

Here the first line represents the clean antiferromagnet, the second is the conventional 
quenched disorder term and the last line contains the static annealed disorder which is 
due to the rare regions. The temperature factor in front of the annealed disorder term 
originates from the Boltzmann factor for the saddle point free energy. The parameter 
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w contains the probability for finding rare regions and the strength of the local order on 
the islands. Since w is non-perturbative in the disorder strength, the theory contains 
effects beyond the conventional perturbative approach. 

This effective translationally invariant theory can now be analyzed by standard 
renormalization group methods, ft turns out |84| ] that the new term in the Landau- 
Ginzburg- Wilson functional ( p9|) destabilizes the critical fixed point found within the 
conventional theory |86[| . No new fixed point is found (at one loop order of the per- 
turbation theory), fnstead, the system displays runaway flow to large disorder values 
in the entire physical parameter space. Within the renormalization group approach it 
is not possible to determine the ultimate fate of the transition. The runaway flow can 
be interpreted either as a complete destruction of the antiferromagnetic long-range 
order in favor of a random singlet phase |8^,pq] or the existence of a non-conventional 
critical point (e.g., with activated scaling). 

While rare regions destroy the conventional critical point in itinerant quantum 
antiferromagnets they do not influence the quantum phase transition of itinerant fer- 
romagnets 184]. The reason is the effective long-ranged interaction between the order 
parameter fluctuations discussed in Sec. 0. It suppresses all fluctuations including 
those generated b y th e rare regions. Therefore the conventional critical behavior dis- 



cussed in Subsec. 3^ will not be changed by the rare regions. 



5 Metal-insulator transitions of disordered interacting electrons 

5.1 Localization and interactions 

Metal-insulator transitions are a particularly fascinating and only incompletely un- 
derstood class of quantum phase transitions. Conceptually, one distinguishes between 
Anderson transitions in models of noninteracting electrons, and Mott-Hubbard tran- 
sitions of clean, interacting electrons. At the former, the electronic charge diffusivity 
D is driven to zero by quenched, or frozen-in, disorder, while the thermodynamic 
properties do not show critical behavior. At the latter, the thermodynamic density 
susceptibility dn/dfi vanishes due to electron-electron interaction effects. In either 
case, the conductivity cr — {dn/dii)D vanishes at the metal-insulator transition. 

The investigation of the disorder-driven metal-insulator transition has a long his- 
tory. Anderson p9] was the first to realize that introducing quenched disorder into a 
metallic system, e.g., by adding impurity atoms, can change the nature of the electronic 
states from spatially extended to localized. This localization transition of disordered 
non-interacting electrons, the Anderson transition, is comparatively well understood 
(for a review see, e.g., Ref. @). The scaling theory of localization Q predicts that 
in the absence of spin-orbit coupling or magnetic fields all states are localized in one 
and two spatial dimensions for arbitrarily weak disorder. Thus, no true metallic phase 
exists in these dimensions. In contrast, in three dimensions there is a phase transition 
from extended states for weak disorder to localized states for strong disorder. These 
results of the scaling theory are in agreement with large-scale computer simulations of 
non-interacting disordered electrons. 

However, in reality electrons do interact via the Coulomb potential, and the ques- 
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Fig. 4 Resistivity of the 2D electron gas a Si-MOSFET as function of carrier concentration 
and doping. The data clearly indicate the existence of a metal-insulator transition, (from 
Ref. @) 



tion is, how this changes the above conclusions. The conventional approach to the 
problem of disordered interacting electrons is based on a perturbative treatment of 
both disorder and interactions (for reviews see, e.g., Refs. ||92,93|). It leads to a scal- 
ing theory and a related field-theoretic formulation of the problem [Q, which was 
later investigated in great detail within the framework of the renormalization group 
(for a review see Ref. ]1^). One of the main results is that in the absence of external 
symmetry-breaking (spin-orbit coupling or magnetic impurities, or a magnetic field) 
a phase transition between a normal metal and an insulator only exists in dimensions 
larger than two, as was the case for non-interacting electrons. In two dimensions 
the results of this approach are inconclusive since the renormalization group displays 
runaway flow to zero disorder but infinite interactions. Furthermore, it has not been 
investigated so far, whether effects of rare regions analogous to those discussed in Sec. ^ 
for magnetic transitions would change the above conclusions about the metal-insulator 
transition. 

Experimental work on the disorder-driven metal-insulator transition (mostly on 
doped semiconductors) carried out before 1994 essentially confirmed the existence of 
a transition in three dimensions while no transition was found in two-dimensional 
Therefore it came as a surprise when experiments on Si-MOSFETs 



96| 



systems. Therefore it came as a surprise when experiments on Si-MOSFETs |95 
revealed indications of a true metal-insulator transition in two dimensions (see Fig. 
0). While these results were first viewed with considerable skepticism they were soon 
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confirmed |97| and later also found in various other materials [Q .[| It soon became clear 
that the main difference between the new experiments and those carried out earlier 
was that the electron (or hole) density is very low. Therefore, the Coulomb interaction 
is particularly strong compared to the Fermi energy. For example, in the Si-MOSFETs 
the typical electron density is lO^^cm"^ leading to a typical Coulomb energy of about 
10 meV while the Fermi energy is only about 0.5 meV. Therefore interaction effects are 
a likely reason for this new metal-insulator transition in two dimensions. A complete 
understanding has, however, not yet been obtained. Different explanations have been 



suggested based on the perturbative renormalization group [ 100 , 101 1 , non-perturbative 
effects 1 103, 103 1, or the transition actuall y bei ng a superconductor-insulator transition 



rather than a metal- insulator transition | 104 ]. In addition to these interaction based 



explanations a number of more conventional suggestions have been made, among them 
the presence of t emp erature-dependent disorder as provided by th e filling and emptying 



of charge traps | 105 | and temperature-dependent screening | 106 ]. 

This is not the place to discuss all these developments in detail or even to review the 
vast field of metal-insulator transitions. Instead, we will concentrate on a few aspects 
of the metal-insulator transition in the presence of both disorder and interactions. 

5.2 Rare regions, local moments, and annealed disorder, a new mechanism for the 
metal-insulator transition 

In this subsection we discuss how rare regions analogous to those studied in Sec. influ- 
ence the metal-insulator transition. Let us consider an electron system in the presence 
of both interactions and (nonmagnetic) quenched disorder. Due to the disorder there 
will be rare spatial regions where the exchange interaction is greatly enhanced. In 
these regions the system will display local magnetic order. Physically, these regions 
correspon d to local magnetic moments. There is much experimental evid ence for local 
moments 



107 1, and their formation has been studied theoretically |108] 



Belitz, Kirkpatrick and Vojta |109] have developed an approach which includes the 
effects of the local moments into a transport theory, starting from a field-theoretical 
description of disordered interacting electrons. As in Sec. the general idea is to 
avoid the disorder average at the beginning of the calculation but to work with a 
fixed disorder configuration. This leads to the appearance of spatially inhomogeneous 
saddle points of the field theory. In particular, there will be saddle points which have 
non-zero magnetization in some rare spatial regions. Analogous to Sec. ^ summing 
over the manifold of degenerate saddle points leads to the appearance of annealed 
magnetic disorder in addition to the underlying (nonmagnetic) quenched disorder. 
Let us emphasize that this annealed magnetic disorder is generically self-generated by 
the system. 

In Ref. |109{ the annealed magnetic disorder was then incorporated into the sigma- 
model description of the metal-insulator transition. To simplify the problem, a model 
of non-interacting electrons with annealed magnetic disorder was considered. This cor- 



Note, however, that very recent experimental data [M on 2d GaAs hole systems indicate that 
the seeming metallic phase is a finite temperature phenomenon. For sufficiently low temperatures 
the old results of scaling theory remain valid, and thus there may be no true metallic phase in two 
dimensions at least in this material. 
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responds to neglecting all interaction effects beyond the formation of local moments. 
The resulting non-linear sigma model can be analyzed using the standard renormal- 
ization group methods |0|. It turns out that the annealed magnetic disorder leads 
to a new mechanism and a new universality class for the metal-insulator transition 
which is different from the conventional localization transition. Note that the effects 
of annealed magnetic disorder are also very different from the case of quenched mag- 
netic impurities. For the simplified model and neglecting the Cooper channel we find 
that the diffusion coefficient is not rcnormalizcd at one loop order while the thermo- 
dynamic density susceptibility dn/dji is driven to zero. Thus, the transition resembles 
a Mott-Hubbard transition rather than an Anderson transition. 

The lower critical dimension for this new transition is two, as it is for the con- 
ventional localization transition. In two dimensions the system is insulating for any 
disorder. Therefore, the local moments alone do not provide an explanation for the 
metal-insulator transition in the two-dimensional electron system in Si-MOSFETs and 
other materials discussed in the last subsection. Clearly, it would be interesting to 



study generalizations of the model studied in Ref. |109| which include the Cooper 



channel and interactions beyond the formation of local moments. 

5.3 Numerical simulation of disordered interacting electrons 

The remaining part of Section g is devoted to the numerical work on interacting 
electrons in the presen ce of que nched disorder. The model investigated is the quantum 



Coulomb glass model [110- I12| , a generalization of the classical Coulomb glass model 
[ 113 , 114 which was used to study disordered insulators. The quantum Coulomb glass 
is defined on a hypercubic lattice of L'^ sites occupied hy N = K L"^ spinless electrons 
{0 < K < 1). To ensure charge neutrality each lattice site carries a compensating 
positive charge of Ke. The Hamiltonian is given by 

H = -tj^icjc, + c]ci) + J2 ^^n^ + ^ ^("^ " ^^^'^^ " ^^^^^^ (^0) 

(ij) i i#J 

where c\ and Ci are the electron creation and annihilation operators at site i, respec- 
tively, and {ij) denotes all pairs of nearest neighbor sites, t gives the strength of the 
hopping term and Ui is the occupation number of site i. For a correct description 
of the insulating phase the Coulomb interaction between the electrons is kept long- 
ranged, Uij — V JTij , since screening breaks down in the insulator (the distance rij is 
measured in units of the lattice constant). The random potential values ipi are chosen 
independently from a box distribution of width 211^ and zero average. Two important 
limiting cases of the quantum Coulomb glass are the Anderson model of localization 
(for Uij = 0) and the classical Coulomb glass (for t — 0). 

For two reasons the numerical simulation of disordered quantum many-particle 
systems is one of the most complicated problems in computational condensed matter 
physics. First, the dimension of the Hilbert space to be considered grows exponentially 
with the system size. Second, the presence of quenched disorder requires the simula- 
tion of many samples with different disorder configurations in order to obtain averages 
or distribution functions of physical quantities. In the case of disordered interacting 
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electrons the problem is even more challenging due to the long-range character of the 
Coulomb interaction which has to be retained, at least for a correct description of 
the insulating phase. Here we discuss the results of two different numerical meth- 
ods to tackle the problem. First, the Coulomb interaction is decoupled by means of 
a Hartree-Fock approximation and numerically solved the remaining self-consistent 
disordered single-particle problem. This method permits comparatively large system 
sizes of more than 10^ sites. The results of this approach are summarized in Sec. ^.4| 
together with those of exact diagonalization studies we performed to check the quality 
of the Hartree-Fock approximation. Since the Hartree-Fock method turned out to be 
a rather poor approximation for the calculation of transport properties an efficient 
method to calculate the low-energy properties of disordered quantum many-particle 
systems with high accuracy has been developed. This method, the Hartree-Fock b ased 



diagonalization, and the results we have obtained this way are summarized in Sec. 5.5 



5.4 Hartree-Fock approximation 

The Hartree-Fock approximation consists in decoupling the Coulomb interaction by 
replacing operators by their expectation values: 

Hhf = - t^iclcj + c]ci) + ^((p, - y))n, 

+ ^ n, [/„(n, -K)-Y, c\cj U,,{c]c,), (41) 

where the first two terms contain the single-particle part of the Hamiltonian, the 
third is the Hartree energy and the fourth term contains the exchange interaction. 
(...) represents the expectation value with respect to the Hartree-Fock ground state 
which has to be determined self-consistently. In this way the many-particle problem 
is reduced to a self-consistent disordered single-particle problem which we solve by 
means of numerically exact diagonalization. 

This method was applied to the three-dimensional quantum Coulomb glass model 



[112|. It was found that the interaction induces a depletion of the single-particle 



density of states in the vicinity of the Fermi energy. For small hopping strength t the 



depletion takes the form of a Coulomb gap [114, 115| known from the classical {t — 0) 
limit. With increasing hopping strength there is a crossover from the nearly parabolic 
Coulomb gap to a square root singularity characteristic of the Coulomb anomaly |116[ 
in the metallic limit. The depletion of the density of states at the Fermi energy has 
drastic consequences for the localization properties of the electronic states. Since 
the degree of localization is essentially determined by the ratio between the hopping 
amplitude and the level spacing, a reduced density of states directly leads to stronger 
localization. Specifically, we calculated the inverse participation number 

^.-'=EKJ»I' (42) 

i 

of a single-particle state \i>) and compared the cases of non- interacting and interacting 
electrons. In the presence of interactions we found a pronounced maximum at the 
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Fermi energy with values above that of non- interacting electrons. Thus, within the 
Hartree-Fock approximation electron-electron interactions lead to enhanced localiza- 
tion. 

In order to precisely determine how the location of the metal-insulator transition 
changes as a result of this effect, we used the fact that the spectral statistics on the 
insulating side is equivalent to that of a Poisson ensemble of random matrices while 
the spectral statistics on the metallic side is that of a Gaussian orthogonal ensemble 



of random matrices [117, IIS]. By analyzing the statistics of the Hartree-Fock levels 
as a function of the hopping amplitude t we could thus determine the location of the 
metal-insulator transition in the phase diagram. In all cases the transition to a metal 
requires a larger hopping strength for interacting electrons than for non-interacting. 

Since the Hartree-Fock method is an uncontrolled approximation it is highly de- 
sirable to compare these results to those of exact many-particle calculations. Now 
the question arises, what quantities are particular suitable for such a comparison. In 
principle, one should compare the values of physical observables like the conductivity 



which is given by the Kubo-Greenwood formula |121, 122] 



0_2 

ReG--(c.) = y2mf\iy)\H{uj + Eo - E,) (43) 

where lo denotes the frequency, j^ is the x component of the current operator and 
\v) denotes the eigenstates of the Hamiltonian. However, the Kubo-Greenwood con- 
ductivity is a complicated quantity from a numerical point of view since it involves 
a nontrivial extrapolation to zero frequency. Therefore, it is desirable to apply the 
simpler localization criteria known from non-interacting systems. This leads, however, 
to an additional complication. The usual criteria which are based on the (single- 
particle) participation number or the statistics of (single-particle) energy levels are 
not even defined for many-particle states. While this does not create any problems at 
the Hartree-Fock level (the Hartree-Fock states are effective single-particle states) it 
implies that the criteria cannot directly be applied to true many-particle states. There 



have been several suggestions (for a discussion see Ref. [119|) how to generalize the 
participation number to a many-particle state. It turns out that a unique procedure 
does not exist, and different generalizations have different merits. Since localization 
is connected with the probability of a particle to return to its starting position, a 
promising candidate to replace the inverse participation number is 

^P (e) = 4t ^ y lim - G,, (e + z(5) G,, {e-ii) . (44) 



3 



Here g{e) is the single-particle density of states and Gjj is the diagonal element of 
the single-particle Greens function. This quantity measures the probability of a single 
particle excitation created at site j to return to this site in infinite time. However, 
in addition to the desired localization information Rp also contains information about 
the decay of the quasi-particles. To disentangle localization and decay information 
we note that the quasi-particle weight zq (which measures the decay only) can be 
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expressed as 

11 A 

^0 (e) = ^-v T7 I] lim - Grj [e + i5) G,, {e ~ i5) . (45) 

The desired localization information is thus contained in the normalized return prob- 
ability P~^ = Rp/zQ which is the natural generalization of the inverse participation 
number ( p2| ) to interacting systems. For non-interacting electrons it reduces to the 
inverse participation number proper. 



Using these criteria the Hartree-Fock results were compared |12C] to those of nu- 
merically exact diagonalizations of the many-particle Hamiltonian. Because the exact 
calculations require the diagonalization of a matrix whose dimension equals the size 
of the Hilbert space this comparison was restricted to systems with not more than 
16 sites. While it was found that the Hartree-Fock approximation yields reasonable 
results for static quantities like the single-particle density of states, it does a very poor 
job for localization properties and for time correlation functions such as the conductiv- 
ity. Thus, the Hartree-Fock results for the influence of the electron-electron interaction 
on transport quantities are highly unreliable. For this reason it was necessary to de- 
velop a new numerical method, the Hartree-Fock based diagonalization which will be 
discussed in the next subsection. 

5.5 Hartree-Fock based diagonalization 

Since numerically exact diagonalizations of the full many-particle Hamiltonian are 
restricted to very small system sizes and, as we have seen, the Hartree-Fock results 
for transport properties are unreliable, a different method which gives exact results 
or at least provides a controlled approximation for comparatively large systems is 



highly desirable. Such a method is the Hartree-Fock based diagonalization |122]. It 



is related to the configuration interaction (CI) approach used in quantum chemistry 



(see, e.g., Ref. [124 ). The basic idea is to diagonalize the many-particle Hamiltonian 
not in a real-space or momentum basis but rather in an energetically truncated basis 
of Hartree-Fock states. The Hartree-Fock states are comparatively close in character 
to the exact eigenstates in the entire parameter space. Therefore it is sufficient to 
keep only a very small fraction of the Hilbert space (e.g., 2000 out of 9 x 10^ states 
for 18 electrons on 36 lattice sites) to obtain low-energy quantities with an accuracy 
comparable to that of exact diagonalization. A schematic of the Hartree-Fock based 
diagonalization method is shown in Fig. pi So far we have carried out test calculations 
for systems with up to 64 lattice sites (32 electrons) and production runs for up to 36 
lattice sites (18 electrons). 

As a first application of the Hartree-Fock based diagonalization method the quan- 



tum Coulomb glass was studied in two dimensions [125|. In order to determine the 
influence of the Coulomb interaction on transport the Kubo-Greenwood conductance 
( [43| ) was calculated. It turned out that the influence of electron-electron interactions 
on the d.c. conductance is opposite in the weakly and strongly disordered regimes. The 
conductance of strongly disordered electrons is considerably enhanced by a weak in- 
teractions. With increasing kinetic energy the relative enhancement decreases as does 
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Fig. 5 Structure of the Hartree-Fock (HF) based diagonaUzation method. 



the interaction range where the enhancement occurs. The conductance of weakly dis- 
ordered electrons is reduced even by weak interactions. In contrast, sufficiently strong 
interactions always reduce the conductance, and the system approaches a Wigner 
crystal or Wigner glass state. These results are summarized in Fig. |^. 

The qualitative difference we found between the weakly and strongly disordered 
regimes can be explained in terms of two competing effects of the interactions. The 
dominant effect in the case of strong disorder is the destruction of phase coherence 
for the single- electron motion. It can be understood from the following simplistic ar- 
gument: Electron-electron collisions are inelastic and thus phase-breaking from the 
point of view of a single electron. Therefore collisions destroy the quantum interfer- 
ence necessary for Anderson localization^ It is clear that this mechanism to increase 
the conductance is particularly effective if the localization length is very small to begin 
with. In the opposite case, i.e., for weak disorder, the dominant effect of the interac- 
tions is an increasing pinning of the electrons by the repulsive forces which for strong 
interactions eventually leads to the formation of a Wigner glass or crystal. This simple 
arguments also explain, why the Hartree-Fock approximation did not reproduce the 
enhanced conductance in the localized regime: The Hartree-Fock potential is static 
and therefore it does not break the phase of the electron - in contrast to real electron- 
electron collisions. 

Analogous investigations were performed for quantum Coulomb glasses in three 



[126| and one |127| spatial dimensions. It turned out that the qualitative disorder and 
interaction dependencies of the conductance are identical to those in two dimensions. 
In particular, it was always found that weak interactions induce delocalization in 
the strongly disordered regime. The degree of delocalization seems to be determined 
essentially by the ratio between disorder strength W and the hopping band width zt 



s Of course, this argument is an oversimplification: Taken at face value, it appears to always 
predict diffusion and thus metallic behavior. In order to obtain a more detailed prediction one has 
to analyze the phase space for electron-electron collisions in the localized regime. Moreover, what 
is neglected in the above argument is that the energy and phase of the many-particle state are not 
changed by a collision. 
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Fig. 6 d.c. conductance of a 2D system of 5 x 5 sites. The disorder strength is W = 1. The 



results represent logarithmic averages over 400 disorder configurations, (from Ref. [125]) 



where z is the number of nearest neighbor sites. In addition to the results obtained 
by the HFD method, other numerical simulations also show a delocalizing influence 
of electron-electron interactions in one [128| and two [110, 111, 129, 130| dimensions. 
The numerical results for the influence of weak interactions on the Kubo-Greenwood 
conductance can be compared to the results from the perturbative renormalization 



group for spinless fermions |131| which predicts that repulsive interactions always 
reduce the transport of disordered spinless fermions. Thus, the numerical results 
agree with those from perturbation theory in the case of weak disorder, while in the 
strongly disordered regime the perturbative results are qualitatively incorrect. 

In addition to the Kubo-Greenwood conductance the above mentioned return prob- 
ability of single-particle excitations was also calculated. The return probability at the 
Fermi energy displays a similar behavior as the Kubo-Greenwood conductance: For 
weak disorder the return probability increases with the interaction strength. The delo- 
calizing tendency in the strongly localized regime is less pronounced, probably because 
the Coulomb gap in the single-particle density of states counteracts the delocalization. 

In summary, the results obtained so far clearly show that interactions can have a 
delocalizing tendency in certain parameter regions. However, it is not clear whether 
this is sufficient to explain the metal-insulator transition found in two dimensions. 
More detailed investigations are presently under way to flnd the flnite-size scaling 
behavior of the conductance. Moreover, we are currently generalizing the simulations 
to include spin degrees of freedom which have been shown to be important in the 
experiments on Si-MOSFETs |l3|,|l3|. 
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6 Summary and outlook 

In this review we have discussed quantum (or zero temperature) phase transitions in 
electronic systems. After a pedagogical introduction on the similarities and differ- 
ences between classical thermal and quantum phase transitions we have considered 
three specific examples in some detail: quantum phase transitions in the spherical 
model, magnetic quantum phase transitions of itinerant electrons, and some aspects 
of the disorder-driven metal-insulator transition. In this final chapter we want to sum- 
marize the results from a common perspective and discuss some of the remaining open 
questions as well as future research directions. 

The theoretical description of a particular phase transition occurring in nature 
usually starts with the identification of the relevant variables, the most important one 
being the order parameter. To proceed further, analytical investigations often follow 
the Landau-Ginzburg- Wilson philosophy, i.e., all degrees of freedom other than the 
order parameter fluctuations are integrated out, resulting in an effective theory in 
terms of the order parameter only. If a rigorous analytical derivation is complicated, 
the Landau-Ginzburg- Wilson theory is sometimes guessed based on general symme- 
try considerations. (This is true in particular when constructing toy models like the 
quantum spherical model.) However, our work on the ferromagnetic quantum phase 
transition has shown that great care has to be taken in such an approach. At zero 
temperature an electronic system usually contains many non-critical soft modes in 
addition to the critical order parameter fluctuations. If these soft modes couple to 
the order parameter integrating them out leads to a singular behavior in the Landau- 
Ginzburg- Wilson theory. This mechanism is not restricted to zero temperature. Any 
soft mode coupling to the order parameter can -- when integrated out - produce sin- 
gularities in the resulting Landau-Ginzburg- Wilson theory. However, in electronic 
systems the number of soft modes at zero temperature is much higher than at finite 
temperatures. 

In the case of itinerant ferromagnets we have shown that these singularities result 
either in unusual non-mean-field scaling behavior at the quantum phase transition or 
in the quantum phase transition being of first order (as is the case experimentally 
for the transition in MnSi). Within the singular Landau-Ginzburg- Wilson theory it 
is very hard to find out what scenario is realized for what microscopic parameters 
since explicit calculations arc next to impossible. A much more promising though 
technically more challenging approach consists in not integrating out all degrees of 
freedom other than the order parameter. Instead, the effective field theory should 
treat all soft modes in the system on the same footing. We will explore the merits of 
such an approach in the future. 

Another potential problem of the standard approach is caused by the fact that 
it is based on perturbation theory: First, integrating out the microscopic degrees of 
freedom in the derivation of the Landau-Ginzburg- Wilson theory can usually be done 
only perturbatively Second, the analysis of the effective theory is done by means of 
the perturbative renormalization group. Now, if the free energy landscape in configu- 
ration space contains many, nearly degenerate minima separated by large barriers, as 
for instance in the presence of random interactions or random fields, a straight-forward 
perturbative approach may miss part of the physics. We have shown this explicitly for 
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the antiferromagnetic quantum phase transition of weakly disordered itinerant elec- 
trons. By approximately taking into account the non-perturbative degrees of freedom 
connected with the many nearly degenerate free energy minima we have derived a cor- 
rection to the effective action of the standard perturbative approach. This correction 
takes the form of static annealed disorder. Including the new term into the renor- 
malization group destabilizes the conventional critical behavior. However, within our 
theory the ultimate fate of the transition could not be determined. Further work will 
be necessary to decide between the different possibilities, viz., a complete destruction 
of the antiferromagnetic long-range order or an unconventional transition that may be 
characterized by an infinite-disorder fixed point and activated scaling. 

It should be emphasized here that we have also shown that including the non- 
perturbative degrees of freedom connected with the many-valley structure of the free 
energy in the presence of disorder does not always destroy the conventional critical 
behavior. For itinerant quantum ferromagncts, for instance, the long-range inter- 
action between the order parameter fluctuations stabilizes the conventional critical 
fixed point. In principle, the method we used to include rare regions into magnetic 
quantum phase transitions of disordered itinerant electrons can be applied to any 
quantum phase transition in the presence of quenched disorder. In the case of the 
disorder-driven metal-insulator transition we found that magnetic rare regions (which 
physically correspond to local magnetic moments) can lead to a new mechanism and 
a new universality class for the transition. In the simplified model investigated so far 
(which neglects all interaction effects beyond local moments) no transition was found 
in two dimensions. An obvious question would be whether this conclusions changes 
for a more realistic model. 

An alternative way to go beyond perturbation theory are numerical simulations. 
We have chosen this approach for an investigation of the metal-insulator transition 
in the presence of both interactions and quenched disorder. We have developed an 
efficient numerical method, the Hartree-Fock based diagonalization method which is 
related to the configuration interaction approach used in quantum chemistry. It allows 
us to calculate with high accuracy the low-energy properties of disordered interacting 
electrons for moderately large systems. So far, we have mostly obtained qualitative 
results concerning the influence of the electron-electron interactions on the transport 
properties: For systems of spinless fermions the interactions reduce the d.c. conduc- 
tance for weak disorder (in agreement with perturbative results) while they enhance 
the d.c. conductance for strong disorder. In future work we plan to extend this work 
to a systematic finite-size scaling study which will lead to quantitative results about 
the metal-insulator transition. Moreover, we are working on including the spin degrees 
of freedom into our numerical scheme. This is of particular importance since the spin 
degrees of freedom have proven to be crucial for the recently discovered metal-insulator 
transition in two dimensions. 
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